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1. INTRODUCTION 

In state control systems, the observers of dynamic systems that are subjected to unknown inputs (Uis) 
are fundamental. In many industrial processes, it is impractical to measure all of the input signals of the system. 
In some cases, the uncertainties of particular system parameters can be described as UIs for the purpose of 
constructing a robust observer. An observer capable of estimating the states in the presence of UIs is required 
for that type of system to develop an appropriate control. Such an issue occurs in systems that are susceptible 
to disturbance or have inaccessible inputs, and it can present in a number of applications like fault detection 
and parameter identification. The main idea behind employing unknown observers for fault detection is to 
generate residuals from the difference between the real and estimated system outputs using an observer. 

The objective of this research is to develop an unknown input observer (UIO) for nonlinear implicit 
systems represented by a discrete-time Takagi-Sugeno (T-S) structure with premise variables depending on 
unmeasurable states, where both the state and output are influenced by the UI. Note that one of the many 
benefits of the T-S fuzzy method is that T-S models provide an excellent technique to describe complicated 
nonlinear systems through the interpolation of many linear models using a membership function. Once the T- 
S models are obtained, the usual linear system theory may be used to analyze the complex systems. Therefore, 
the T-S models become powerful technical tools for modeling, design observer and control systems, as shown 
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in [1]-[5], and the references contained therein. Otherwise, implicit systems, likewise called descriptor or 
singular systems are governed by a combination of differential equations and algebraic ones and have been 
extensively utilized in dynamic process modeling to represent the behavior of numerous industrial systems, 
e.g. electrical systems, robotics, and chemical processes. Some real utilizations of implicit systems may be 
found in [6]-[8]. 

The fuzzy UIO design problem for explicit or implicit T-S models attracts more and more attention 
and is still largely an open problem and remains an area for additional study in both continuous-time and 
discrete-time cases. Indeed, different advancements in the observer and its deployment to fault detection are 
discussed in the literature. A novel fault detection method for a robot manipulator is provided in [9], which 
utilizes a high-order sliding mode observer to adjust for both uncertainty and faults. A fuzzy observer design 
approach using T-S systems to accomplish the estimation of vehicle state taking into account of the longitudinal 
lateral coupling dynamics and the presence of UIs is developed in [10]. Developing an UIO for T-S systems 
by using an expanded state vector is the aim of the work described in [11]. Bouassem et al. [12] suggested 
another technique based on the separation between static and dynamic relations in the T-S implicit system to 
estimate both the system state and the UIs concomitantly. This approach allowed the estimation of both non- 
measurable states and unknown faults in the actuators and sensors for a class of continuous-time T-S implicit 
models in [13]. Louzimi et al. [14], a design of an UIO is proposed for a class of continuous-time T-S implicit 
systems satisfying the Lipschitz conditions. Similarly, in the case of discrete-time T-S models, various 
publications published for explicit and implicit structures such as [15]-[20]. Essabre et al. [21] suggest a 
technique for constructing an observer for a class of T-S implicit systems subject to UIs utilizing measurable 
premise variables that are dependent on the measurable states, input, output, or external parameters. Yoneyama 
[22], it is confirmed that even if the output of the system is selected as a premise variable and this output is 
directly impacted by perturbations (which seems to be the case in majority of realistic conditions), the T-S 
system acquired does not represent an exactly nonlinear system. For this reason, using the system state as a 
premise variable makes it possible to depict a broader class of nonlinear systems. 

Motivated by the above discussions, the purpose of our study is not to compare our technique to 
previous ones, but rather to generalize our results from the case of T-S implicit models with measurable premise 
variables [21] to the case of T-S implicit models with unmeasurable premise variables. Using the singular value 
decomposition technique, the fundamental contribution of this research is to develop a UI observer design for 
a class of discrete-time Takagi-Sugeno implicit systems (DTSISs) models with unmeasurable premise 
variables, allowing for the estimate of UIs and unknown states at the same time. In comparison to T-S models 
that employ just measurable premise variables, the T-S model using unmeasurable premise variables is an 
important structure because it can represent a wider variety of nonlinear systems and can describe precisely the 
behaviour of the general nonlinear system. The convergence of the state estimation error is established by using 
the Lyapunov theory, and the stability criteria are expressed as linear matrix inequalities (LMIs). The suggested 
result in this research is presented in an explicit form. 

The paper is organized as follows: the mathematics of the DTSISs subject to UI to be researched is 
explained in section 2. Section 3 introduces the major result of this work. Finally, this theoretical result is 
applied to a model of the rolling disc to show the excellent performance of the obtained result. 


2. DESCRIPTION OF THE CONSIDERED MODEL 
The model of the DTSIS that we study in this article is of the following form: 


poe = Wh Pilk) (Aire + Biur + Cidg) (1) 
Yk = Dxx + Edk 
where x; € R”: the states vector, Ug € R™: the control input, ye R?: the measured output vector, dge R”: the 
unknown input Me R”*”, The matrices such that rank(M)<n, A; € R"“", B; €e R®™, C; €e R®*T, D e RP*" and 
E e RP*" are real known constant matrices. 
The premise variable is included in the vector x. y;(&,) are the weighting functions that guarantee 
the transition between the contributions of each sub-model: 


ae = AjX, + Byu, + Cidk 2 
VR = Dxx + Ed, ( ) 

they hold the convex su property: 
we GE) = 1; 05 9;(&,) $1; (3) 
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i= 1,...,q 


we consider in this study that the weighting functions p; (x) are dependent on unmeasurable premise variables 
(system state). 


2.1. Assumption 1 
Assume d is taken as a constant in each interval, i.e.: 


dua = dy; k € [T, T2]; Y T T, E€ R* (4) 


to concomitanly estimate the UIs and system state for system (1), we use the augmented system structure, 
defined as follows: 


ae = YE, gin) (Qing + Siuk) (5) 
Ve = Fnk 


where: 


mt) . 

ea c) (6) 
Si = (9 

F=(D E) 


Before we get to the main point, consider the subsequent assumption. 


2.2. Assumption 2 
Suppose that: i) H1) (H, Q;) are regular, i.e. det(zH — Q;) # 0 Vz E C, ii) H2) All sub models (2) 
are impulse observable and detectable, and iii) H3) rank ($) =0 =Nn+r. 


Taking into account the hypothesis H3), there exists a non-singular matrix C A such that: 
3 J4 


LH+}F =! 
ai +a = 0 (9 


where J, € R°*?, J} € R°*?, J, € RP“? and J, € R?*? are constant matrices determined from the singular value 
decomposition of (5). 


3. MAIN RESULT: FUZZY UNKOWN INPUTS OBSERVER DESIGN MAI 

The goal of this paragraph is to provide a novel observer design approach to estimate the UIs and the 
unmeasurable states of DTSIS (1). To do this, the DTSIS (1) is converted into the equivalent system (5). Thus, 
the suggested fuzzy observer has the following structure: 


aa = EL pilk) (Nilik + Lite + LziYk + Giu) (8) 
Ae = Uk +J2Yk + Kaye 


Where f},: the estimated augmented state vector; ĉẹ: the estimate of the premise vector, the gains matrices N;, 


Lii, Lz; and G;must be calculated such that 4, converges to ng exponentially. Let us develop the augmented 
state estimation error in order to set the conditions for the exponential convergence of the observer (8): 


Ek = Nk — Tk (9) 
by substituting (5), (7) and (8) into (9) we obtain: 


& = (Ki, + KJ3Yk)HNk — Hk (10) 
Design of unknown input observer for discrete-time Takagi Sugeno fuzzy implicit ... (Mohamed Essabre) 


62 o ISSN: 2302-9285 
from (5) and (8) it follows that the dynamics of this error is: 
Est = Lier Pil) Ua + KJ3) (Qik + Site) — Dey Pilk) (Mite + Lite + Laie + Gi.) UD 


using (10), (11) can be written as: 


aa ` Pi (EU + KJ3) (Qiy + Sitix) + ` pi (Êk) Nir - J PENNO ay 


+ KJ3)H + LiF + LziF)nk + G; jUx) 
provided the matrices N;, Lii, L2; and G; and K satisfy: 
N;Qı + KJ3)H + LiF + LaF = Uy + KJ3)Q; (13) 
Gi = (Q1 + KJ3)S; (14) 


then, using (7) and (13), we get: 


N; = Gi + KJ3)Q;i — LaiF + (Ni U2 + KJ4)— La) F (15) 
take: 

Ly; = N;U2 + KJ4) (16) 
then: 

N; = G1 + KJ3)Qi — LiF (17) 


it follows system (12) is equivalent to: 


Exar = Diaa Pilé) Nier + LL, (Gi) — pila + KJ) (Qing + Situr) (18) 


note that: 


YL, iE.) — Pil = ye PilEx) P; Ex) (Qi = 05) (19) 


Yoto- gE NS = y PENEN: - S) 


i j=1 
the (18) is then transformed as (20): 
AGAMIA + K]J3)(AQijnx + AS; Ux) (20) 
where AQ; = Q; — Qj and AS;; = S; — S;. By multiplying by x ij= 1(Pi(E), (20) can be simplified to the (21): 
Ek+1 = Ei j1 PilEx) Pi (Ex) (N jEr + Dine + Tjur) (21) 


where 


®ij = Gh + KJ3)AQij 
ij €(1,..,q) 


considering & = (ef ni)", we get: 
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H êk = Xo PiE) 9; (E) (QijE + Siju) (23) 
where 
y I 0 
= (5 A 
_ IN, D; 
Be) (24) 


To introduce the convergence conditions of the suggested observer (8), we form the following 
theorem. Theorem: according to hypotheses 1 and 2, if for DTSIS (1) there exist symmetric matrices P,, P2, Q, 
and the matrices W;, j € (1,...,q) satisfying the LMI (25), then the observer gains that guarantee the 
convergence to zero of the estimation error can be determined: 


T T 
Miq M21 M31 Magy 
T T 
m m m m 
21 22 32 B < 0 (25) 
M31 M32 M33 Maz 
Ma, Maz M43 M44 


where 


mıı = -P 
M22 = —Q} PzQi — H" P,H 
M33 = SÍ P25; 


Ma, = —Py 
M2, = 0 
M3, = 0 2s) 


Maz = PyJ1Q; + QJ3Q;-W;F 
Maz = PJ AQij + QJ34Qij 
M32 = Si P,Q; 

M43 = 0 


with J; , Jz , Js, and J, are such that (7) is satisfied. 
The observer gains N,, L1;, and G; are given such that (14), (16), and (17) are satisfied. K and L2; are 
obtained by: 


K = Pe 

a Š y ; (27) 
Proof: the following quadratic Lyapunov function is used to demonstrate the theorem: 

V; = &' A PHé,, ATPH 20, P= P™>0 (28) 
with: 

P= (es 3 ) (29) 
its variation AV = V,,, — Vg along the error dynamics (23) is calculated by: 

AV = êk41 HA PH&4, — êk HPH êk (30) 


by using (23), (29) becomes as (31): 
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AV = Dai PilE) 0 (Ex) Qije + Siju)" P(Qijé + Sjur) — êk A PAE, (31) 


multiplying by Èi j= Pip; (Êk), (30) can be reduced to the (32): 


ar é 
AV = Fons GOE) e uh] My || (32) 
where: 
2 a - AT PA ars) es 
SUP Qi SHPSij 
Lj E€ (l4) 
thus, the convergence of the estimation error is ensured if the following condition is guaranteed: 
OX; <0 (34) 
from (6), and (24), (34) can be written as: 
911 942 43 
021 O22 823} <0 (35) 
931 032 O33 
ij,€ (,..,q) 


011 = NF P N; -P 

62, = pi PN; 

031 = Ij PaM; 

642 = Nj P, Gi; 

022 = PiP Pi; + Qi P,Q; me H™P,H (36) 
032 = Tj Py iy + ST P2Q; 

013 = Nj Pl 

023 = PpP Ij + Q;P2S; 

033 = Dj Ply + Si PzSi 


based on the schur complement [23], (35) is equivalent to the following conditions: 


=e 0 0 NIP, 
0 QIP Qi— HP H QPS; PiP Sô (37) 
0 Si P2Qi SPP2S; ij Py 

PAN P®y Ply -P 

ij,€ (1, ...,q) 


then, from (17), and (22), and use of the change of variable: 


Q = PK 
fw = PiL2j 2) 


the LMIs conditions (25) are obtained. The proof is completed. 
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4. NUMERICAL APPLICATION TO ROLLING DISC 
To evaluate the suggested observer, we consider the following DTSIS of the rolling process given in 
[24] which is assumed to be affected by an unknown input variable as follows: 


aaa = Vier Pik) (Aire + Biuk + Cd) (39) 
Yk = Dx, + Edk 

where x, = [Xik , =; X4pg |: the state vector, x,,: the position of the centre of the disc, x2: the translational 

velocity of the centre of the disc, x3,: the angular velocity of the disc, x4;,: the contact force, ug and yg: the 


control and the output measurements vectors respectively. dg: the unknown input. The numerical matrices 
values are: 


1 0.01 0 0 1 0.01 0 0 
A, =| 70.025 0.9925 0 0.0003 | 4 _( -0.027 0.9925 0 0.0003 
0 0.01 -0.008 0 p°? 0 0.01  -0.004 0 
—0.025 -0.0075 0 0.0008 —0.025 —0.0075 0 0.0008 
0 1000 0 
1011 
pa Oe c=(0 a4 o), ma[9 19 0) -gahto 
0 AA 00 0 0 0 
—0.0013 00 0 0 0 


The weighting functions are given by: ~,(é,) = 1—12.7551x?, and g,(&,) = 12.7551x2,. Note that to 
apply the suggested UIO (8) for (39), just rewrite the model (39) in its equivalent form (5) as indicated in 


section 2. Therefore, under hypothesis H3) the following matrices J, , J2 , J3 and J, verifying (7) have been 
obtained: 


0.8 02 0 0 0 0.2 —0.2 —0.2 
0.2 08 0 0 0 —0.2 0.2 0.2 
J, =|] -0.2 02 0 0 0j, Jz =} 0.2 0.8 —0.2 
—0.4 -06 0 0 0 04 —0.4 0.6 
0 0 0 0 1 0 0 0 
—0.2828 0.2828 -—0.5774 —0.5164 0 
h= ( 0.2 —0.2 —0.8165 0.3651 o) 
0.2828 —0.2828 0 —0.7746 0 


—0.2 0.2 0.2 


0.2828 —0.2828 —0.2828 
I, = ( ) 
—0.2828 0.2828 0.2828 


the following observer gains result from the solving of LMIs (25) of the Theorem: 


0.0355 0.0355 0 —0.0355 0.0265 
0.0355 0.0354 0 —0.0354 0.0269 
N, = 0 0 0 0 0.0002 
—0.0355 —0.0355 0 0.0355 —0.0269 


0.0258 0.0232  —0.0009 —0.024 0.3628 


0.0355 0.0355 0 0.0355 0.0265 
0.0355 0.0354 0 —0.0354 0.0269 
No = 0 0 0 0 0.0002 
0.0355 —0.0355 0 0.0355 —0.0269 
0.0258 0.0232 —0.0009 —0.024 0.3628 
~1.4729 14729 1.4374 —1.4729 1.4729 1.4374 
—1.5011 15011 1.4657 —1.5011 1.5011 1.4657 
Lı =| -0.0155 0.0155 0.0155 |,Li,=] -0.0155 0.0155 0.0155 
1.5020 —1.5019 —1.4665 1.5020 1.5019 —1.4665 
~22.9698 22.9689 22.9457 ~22.9698 22.9689 22.9457 
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—1.6832 2.0345 1.7173 —1.6832 2.0345 1.7173 
—1.727 2.0782 1.7611 —1.727 2.0782 1.7611 
Lo, = | —0.0242 0.0241 0.0242 |, L22 = | —0.0242 0.0241 0.0242 
1.7285  —2.0797  —1.7626 1.7285  —2.0797  —1.7626 
63.6925 —63.4787 —63.6366 63.6925 —63.4787 —63.6366 
0.0033 0.0033 55.5406 68.3030 —1.4106 
0.0034 0.0034 54.7603 68.8198 —0.5902 
G, = 0 i G, = 0 , K =] —0.2466 0.1520 0.2673 
—0.0034 —0.0034 —55.0423 —68.6124 0.8740 
—0.0796 —0.796 —92.2846 130.4449 40.7553 


The UI, which is applied during the interval [2s, 4s], is shown in Figure 1. 

The simulation results shown in Figures 1-4 illustrate that the presence of the UI causes an undesirable 
and noticeable variation in system states at the moment of t the UI appears. Consistent with the results from 
Figures 1-5, we can see that the suggested observer developed by employing these numerical gain values N;, 
Lii, Lz; and G; permits the estimated states to converge to the actual state as rapidly as possible, restoring 
desired performance and ensuring stability despite the existence of UI. 

To show the advantage and justify the effectiveness of the suggested method, we compare the 
proposed estimator with other advanced methods. For example, the technique described in [25] is founded 
totally on the separation between the differential equations and algebraic ones inside the DTSIS that is 
applicable only for systems that satisfy the criteria of rank specified in the study. Another benefit of our 
technique over others that employ DTSIS with premise variables satisfying Lipschitz criteria (see [14]) is the 
absence of the Lipschitz assumption of the weighting functions, that allowed us to reduce the LMI constraints. 


t{s) 


Figure 1. Unknown input d and its estimate 


ò 
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3 
t{s) t(s) 
Figure 2. State variable x, and its estimate Figure 3. State variable x, and its estimate 
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i i 4 4 2 


o 1 2 3 = 5 6 3 4 3 6 
t(s) t{s) 
Figure 4. State variable x3 and its estimate Figure 5. State variable x, and its estimate 
CONCLUSION 


In this research, a new UIO design method for DTSIS with unmeasurable premise variables is 


presented. The approach is based on the singular value decomposition and the utilization of an augmented 
system composed of the state variables of the model and UIs. The suggested solution allowed for the estimation 
of both the system status and the UIs at the same time. The Lyapunov theory is used to investigate the 
convergence of the state estimation error, and the presence of the condition assuring this convergence is 
represented in terms of LMI. Finally, an example of application to a rolling disc is presented to verify the 
performance of the suggested technique in estimating the states and the UI. 
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